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Abstract. We show that the low energy behaviour of quite diverse impurity systems can be described by a 
single renormalized Anderson model, with three parameters, an effective level id, an effective hybridization 
V, and a quasiparticle interaction U. The renormalized parameters are calculated as a function of the bare 
parameters for a number of impurity models, including those with coupling to phonons and a Falikov- 
Kimball interaction term. In the model with a coupling to phonons we determine where the interaction of 
the quasiparticles changes sign as a function of the electron-phonon coupling. In the model with a Falikov- 
Kimball interaction we show that to a good approximation the low energy behaviour corresponds to that 
of a bare Anderson model with a shifted impurity level. 

PACS. 71.10.-w Theories and models of many-electron systems - 71.27.+a Strongly correlated electron 
systems; heavy fermions - 75.20.Hr Local moment in compounds and alloys; Kondo effect, valence fluctu- 
ations, heavy fermions 



1 Introduction 

The low energy behaviour of impurities in a metallic host 
can be calculated using a numerical rcnormalization group 
(NRG) approach in which the higher energy states are 
progressively eliminated |1I2| . The low temperature ther- 
modynamics of the impurity are deduced from the leading 
corrections to the low energy fixed point of the renormal- 
ization group transformation. An alternative approach to 
the calculation of the low energy behaviour is the renor- 
malized perturbation theory (RPT) When this ap- 
proach is applied to the Anderson impurity model, for 
instance, the effective low energy model generated is just 
a renormalized version of the Anderson model; expres- 
sions for the renormalized parameters can be derived in 
terms of the self-energy, its derivative and the irreducible 
four- vertex, all evaluated at zero frequency |3I4| . In this 
paper we will clarify the relation between these two ap- 
proaches, and consider the various possible ways of deter- 
mining the renormalized parameters. We will also general- 
ize the approach to include models which have a coupling 
to phonons and a Falikov-Kimball screening interaction 
between the impurity and conduction electrons. 
The Anderson model has the form, 



matrix element Vk ■ There is in addition a local interaction 
U between the electrons on the impurity site. When U — 
the local level broadens into a resonance, corresponding 
to a localized quasi-bound state, whose width depends on 
the quantity A(ui) — 7r^ fe | Vfc | 2 <5(cj — e^). It is usual to 
consider the case of a wide conduction band with a flat 
density of states where A(u) becomes independent of u 
and can be taken as a constant A. 

In the renormalized model used in the RPT approach 
there is an effective energy level at id, and an effective 
resonance width A. Expressions for these renormalized 
parameters can be derived in terms of the 'bare' parame- 
ters, Ed, A, the local self-energy S a (iv, h) and its frequency 
derivative E' a (oj,h) evaluated at zero frequency w = 0, 
zero magnetic field h — 0, and T = 0. They are given by 



id = z(ea + 27,(0,0)), A = zA, 



(2) 



where z, the wavefunction renormalization factor, is given 
by z = 1/(1— 27^(0,0)). The effective local interaction U 
is expressed in terms of the local four- vertex Jf ^(wi, uj^, u>3, W4) 
evaluated at zero frequency (u>i = L02 = 0J3 = 0J4, = 0), 



Ham = e dd\d a + U n dA n dA + ^ £k,o4.^c k>a + 

(j k,<7 

+ J2(V k dU kt<T + V k *c{ er d (T ). 



i> = z 2 r u (o,o,o,o). 



(3) 



(1) 



The renormalized effective model has the same form 
as Q, but in terms of renormalized parameters, and the 
It describes a localized level e<2 of an impurity, hybridized interaction term is normal-ordered, as it only comes into 
with the conduction electrons of the host metal via the play when two or more excitations are created relative to 
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the ground state of the interacting system, 



(4) 

where the colon brackets indicate that the expression within 
them must be normal-ordered. This renormalized model 
is similar to that used in earlier phcnomcnological local 
Fermi-liquid theories but here it also includes a quasi- 
particle interaction term. 

In the renormalized perturbation theory the Hamilto- 
nians of the bare and renormalized Anderson models are 
related via Ham = Ham + H c , where H c is the counter- 
term Hamiltonian given by 



which corresponds to the Friedel sum rule [5]. These two 
results correspond to a local Fermi-liquid theory with non- 
interacting quasiparticles and can be deduced from the 
zero order (U = 0) renormalized model. The temperature 
dependence of the impurity contribution to the conductiv- 
ity cr(T) has also been calculated for the symmetric model 
to order T 2 from the renormalized self-energy calculated 
to second order in U. The result is 



H c = ^A x 4<k + A2 



nd,t n d,l- 



(5) 



The renormalized perturbation expansion is in powers of 
the renormalized interaction U, but all the terms in Q 
and jnj are taken into account; the parameters Ai and A2, 
and a rescaling factor A3, are determined by the condition 
that there is no further renormalization of the already 
fully renormalized parameters, e d , A and U arising from 
the expansion (the procedure is clearer in the Lagrangian 
formulation, for details see 0]). It has been shown $ that 
the normal-ordered renormalized Anderson model 0} is 
sufficient for the calculation of thermodynamic properties 
in the low energy, low field, and low temperature regime. 
The counter-term part of the Hamiltonian need only be 
taken into account explicitly for the calculation of correc- 
tions for higher temperatures and magnetic fields jlj . The 
first order expressions for the impurity spin and charge 
susceptibilities, Xs,imp and Xc.imp at T = derived from 
the renormalized model are exact [J] and given by 



ISjimp 



oPimp(0)(l 



Up imp (0)), 



Xc,imp = X/5imp(0)(l - ^p imp (0)), 



(6a) 



(6b) 



where the spin susceptibility is given in units of (,9/ib) 2 
and the charge susceptibility differs by a factor of \ from 
the usual definition, so it is the isospin equivalent of Xs,im P - 
The quasiparticle density of states at the Fermi level pi mp (0) 
is given by 

Pimp(0) = (7) 



e d + A 2 



Exact results for specific heat coefficient 7i mp and the oc- 
cupation of the impurity level rii mPj(T are 



and 



_2tt 2 . 

7imp — g Pimp (.U J, 



2 7T \A 



(8) 



(9) 



u 



0(T 4 ) 



( 10 ) 

When the renormalized parameters are expressed in terms 
of the self-energy and the vertex function these results co- 
incide with the exact expressions derived by Yamada [5] 
from an analysis of perturbation theory to all orders U, 
and the Fermi-liquid results of Nozieres in the Kondo 
limit. This result has also been generalized to include the 
leading non-linear correction to the differential conduc- 
tance through a quantum dot in the Fermi liquid regime 
[TT] . which is of order V e 2 , where voltage V e is the volt- 
age difference across the dot. However, the values of the 
renormalized parameters e d , A and U , have to be deter- 
mined. In the localized or Kondo limit they can be re- 
duced to one single parameter, the Kondo temperature 
Tk- For U 3> \ed\, and e d -€1 0, the electron in the d- 
state at the impurity site will be localized and as a con- 
sequence the impurity charge susceptibility must vanish 
in this limit. From equations, © and J^J, this implies 
U = ttA and e d = 0. If the impurity susceptibility at 
T = is expressed in terms of a Kondo temperature Tk 
defined by Xs,imp = 1/4Tk, then the renormalized param- 
eters can be expressed in terms of a single energy scale 
Tk, so U = it A = 4Tk, and the Wilson or \H ratio, 

r= i+up d (o) = 2 rmj. 

It is clear that the renormalized Anderson model de- 
rived in the RPT must be directly related to the low en- 
ergy effective model obtained in the numerical renormal- 
ization group (NRG) calculations by Wilson [jQ for the 
s-d (Kondo) model, and by Krishnamurthy, Wilkins and 
Wilson (KWW) for the Anderson model. The RPT and 
NRG approaches are, however, rather different. In the Wil- 
son approach the conduction band of the Anderson model 
is replaced by a discrete spectrum of states, which is then 
expressed in the form of a tight-binding chain with the im- 
purity at one end. The Hamiltonian for the discrete model 
for a finite chain with N + 2 sites, including the impurity 
site, is 



rrN 
11 AM 



Un dA n d .i + VJ2( d t c o,a + 4,<A) 



+ J2 A ~ n/2 Zn(cl a C n+ha + 4 +1)ff C„ iCT ), (11) 
n—0,(T 

where A > 1 is the discretization parameter, and £ n is 
given by 



D (l + yl- 1 )(l -A- 11 - 1 ) 

~2 (1 - yl-2n-l)l/2(l _ ^-2^-3)1/2 ' 



(12) 
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and 2D is the width of the conduction band pQ. The dis- 
cretization of the conduction band is logarithmic and such 
that the density of levels increases as the Fermi-level is 
approached, and e = is the limit point of the sequence. 
With the model in this form an iterative diagonalization 
scheme is then set up, starting with the impurity site, and 
an extra site along the chain is added at each iteration 
step. After five or six steps the matrices become too large 
to handle, and after this point the Hibert space is trun- 
cated, with the higher energy states being neglected. A 
fixed number, of the order 500-1500 of the lowest lying 
energy states, is retained at each subsequent step. The 
couplings along the chain fall off as yl~"/ 2 , where n is 
the nth site as A > 1. Lower and lower energy scales are 
reached in the process but the Hamiltonian is rescaled af- 
ter each step by a factor A 1 / 2 , so that the lowest energy 
scale is formally the same after each iteration. 

A renormalization group transformation can be set up 
in which the states and couplings can be compared after 
each step. This renormalization group transformation for 
the s-d (Kondo) and Anderson model has a low energy 
fixed point corresponding to states of a free chain uncou- 
pled from the impurity. In the s-d case this was interpreted 
as a J — > oo fixed point, so that the impurity and first site 
become essentially uncoupled from the chain. As this fixed 
point Hamiltonian corresponds to a free chain it gives no 
finite contribution to the impurity susceptibility and spe- 
cific heat. The finite impurity susceptibility and specific 
heat arise from the leading irrelevant terms of the renor- 
malization group transformation, which correspond to a 
residual interaction and a one-body coupling term at the 
end of this decoupled chain. 

Similar results hold for the Anderson model [2] which 
for the symmetric model can be interpreted as a V — > oo 
fixed point. However, from the Bethe ansatz exact solution 
of the Anderson model |12ll.{| we know that the Ander- 
son model behaves as a Fermi-liquid, whatever the value 
of the interaction term, so the nature of the fixed point 
is independent of U. It seems appropriate from a Fermi- 
liquid point of view, therefore, to base the interpretation 
of the low energy fixed point on the original chain includ- 
ing the impurity, but with renormalized parameters. In 
this approach the fixed point and leading irrelevant terms 
are those of the renormalized Anderson model; this then 
makes a direct connection with the effective model used 
in the renormalized perturbation theory. It also has the 
advantage that it can be interpreted in terms of a quasi- 
particle picture in 1-1 comparison with the original model. 
We can use the NRG technique to calculate the renormal- 
ized Fermi-liquid parameters. A knowledge of how these 
renormalized parameters, such as the quasiparticle inter- 
action U , depend on the bare parameters of the model can 
give us considerable insight into the low energy behaviour 
of impurity systems. The numerical results for the physical 
behaviour of the model at low temperatures, the suscep- 
tibilities, specific heat coefficient etc, are essentially the 
same as those obtained in the earlier NRG calculations of 
KWW the difference here is that they are interpreted 
in terms of a single renormalized model, with parameters 



in 1-1 correspondence with those of the original Ander- 
son model, and the calculational procedure has also been 
considerably simplified. The approach can be extended to 
more general impurity models including orbital degener- 
acy |T3| . 

In the next section of this paper we give a new way 
of analysing the fixed point to determine the renormal- 
ized parameters, which is straight forwardly applicable to 
systems with a non-symmetric and non-constant density 
of states. The approach is then applied to generalizations 
of the model to include an interaction with phonons, and 
Falikov-Kimball terms, for which there are in general no 
exact results. We examine how the quasiparticle interac- 
tion and renormalized parameters vary according to the 
parameter regimes of these more general impurity models. 
In the conclusion we discuss the question of the extent to 
which these additional interactions, such as those due to 
phonons and the Falikov-Kimball terms, may be absorbed 
by modifying the parameters of the bare model, or whether 
to describe fully the behaviour in higher energy regimes, 
these terms have to be included explicitly. 



2 Calculation of Renormalized Parameters 

In this section we will examine the low energy NRG fixed 
point of the Anderson model as a renormalized version of 
the same model and deduce the renormalized parameters 
id, A and U . The starting point of the NRG calculation is 
the discretized form of the model given in equation ljll(l . 
The many-body states are calculated using the iteration 
procedure as outlined in the previous section. For a given 
N, we denote the minimum energy required to add a sin- 
gle electron to the ground state by E P (N), and and the 
minimum energy to create a single hole by Eh(N). If these 
single particle and hole excitations correspond to a renor- 
malized Anderson model, then asymptotically they should 
coincide with those of the free model with an appropriate 
choice of id and A (the term in U plays a role only when 
more than one single particle excitation is created from 
the ground state). This problem is considered in the ap- 
pendix where we define A^-dependent quantities, id(N), 
and A(N) = ttV(N) 2 /2, via the two equations, 



nE p (N)A 



-(N-l)/2 



2A(N) 



2A{N) 



-TrE h (N)A-( N -V/ 2 ir~e d (N) 
2A{N) 



2A(N) 



A^/ 2 g 0Q (E p (N)), 
(13) 

A^/igooi-E^N)). 



The renormalized parameters id and A corresponding to 
the low energy fixed point are given by id = limjv— »oo ^d(N) 
and A = linijv^oo A(N) and are given by equations (|41|l 
and 142(1 in the appendix. 

In figure 1 we plot the quantities A(N) and id(N) 
against TV for a model with U/ttA = 6.0 and ed/^A = 
—4.0. For the impurity site N = we take A(0) = A, and 
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ed{0) to be the average Hartree-Fock value + Und/2, 
where rid is the total occupation value at the impurity site, 
calculated within the Hartree-Fock theory. This Hartree- 
Fock value for €d(0) is considerably shifted from the bare 
value e<i, but clearly fits the trend in values for id(N) for 
small N. If the one-particle excitations can be described by 
an effective non-interacting Anderson model then A(N) 
and £d(N) should be independent of N. This can indeed 
seen to be the case for N > 40 so for this range these 
excitations can be described by a model with renormalized 
parameters A and id- As the bare parameters correspond 
to a model in the Kondo regime one sees that in this case 
Id ~ 0, whereas A/A is small but finite, A being of the 
order of the Kondo temperature Tk- 

Once the renormalized parameters id and V have been 
determined the free quasiparticle Hamiltonian can be di- 
agonalized and written in the form 

(N+2)/2 

A -(N-i)/2 J- (E Ptk (N)pl ia p kta + E h , k (N)hl a h ktir ), 

k=l 

(15) 

are the creation and anni- 
hilation operators for the quasiparticle and quasihole ex- 
citations, and A-( N -V/ 2 E pM {N) and A-^ N -^/ 2 E h>k (N) 
are the corresponding excitation energies relative to the 
ground or vacuum state |0); the scale factor y l~( Ar ~ 1 )/ 2 is 
due to the fact that the energies are calculated for the 
rescaled Hamiltonian, which is such that E Ptk {N) and 
Eh,k(N) for k = 1 are of order 1. For the lowest-lying 
level particle and hole levels, we have E P (N) = E p i(N) 
and E h {N) = E h ^(N). 

From the parameters, id and A, which determine the 
free quasiparticle excitations we can immmediately deduce 
the occupation of the impurity level 7ii mp at T = from 
equation 0, which depends only on the ratio id/ A, the 
quasiparticle density of states from J7J, and the impurity 
specific heat coefficient 7 imp from JHJ). 

To calculate the spin and charge susceptibilities and 
the low temperature dependence of the conductivity we 
need to include a quasiparticle interaction term, which for 
the rescaled model takes the form, 



where p\ p k and h\ , h k 



H U (N) = UA^-^ 2 



d^d^d^d^ : 



(16) 



Asymptotically as TV — * oo the effect of this term on the 
low-lying many-particle excitations tends to zero so it is 
sufficient to calculate the effect of this term to first order 
in U. To this end we need the operator d a expressed in 
terms of the eigenstates of 1)1. 5(1. 



(N+2)/2 
k=l 



■l)pka + M-^)hD- ( 17 ) 



If the lowest two-particle excitation from the ground state 
for the interacting system for a given N has an energy 
Epp(N), then we can calculate U by equating the energy 
difference E pp (N) — 2E p (N) to that calculated using 
asymptotically in the limit N — > oo. For finite N we can 






& a a <B- 




A A A A A A 








,,!,,,,],,, 



& <8 ft O ^ 



ft s 



0.04 
0.03 
0.(12 
0.01 
o 

-0.01 



Fig. 1. Plots of the parameters, A(N)/irA (crosses), 
i d (N)/nA (triangles), U pp {N))/irA (diamonds), U hh (N)/irA 
(stars) and U p h(N)/wA (circles), for the model with bare pa- 
rameters, U/nA — 6.0 and ed/nA = —4.0. 



use this equation to define an W-dependent renormalized 



interaction U pp (N), 

E PP (N)-2E P (N) = U(N)A( N -W 2 \r p ,A-l)\ 2 \r p ,x(- 
where \ip P: i(— 1)| 2 is given by 



l-V 2 (N)A("-i)g> 00 (E p (N)) 



-1)| 2 , 
(18) 



(19) 



where g' 00 (uj) is the derivative of goo (to). 

Alternatively we could consider the same procedure for 
a two hole excitation E^^N) and in a similar way define 
an A^-dependent renormalized interaction Uhh(N), or a 
particle- hole excitation E p h(N) to define a renormalized 
interaction U p h{N). In this latter case, as a positive U 
leads to particle-hole attraction, we use E P (N) + Eh(N) — 
E p h(N) on the left-hand side of equation |JTHJ|. 

If these two particle excitations can be described by 
an effective Anderson model then U pp (N), Uhh(N) and 
U p h{N) should be independent of N and also independent 
of the particle- hole labels. In figure 1 we also plot U PP (N), 
Uhh(N) and U p h(N) as a function of A^ for the parameters 
used earlier, U/irA = 6 and ed/nA — —4, with values at 
N ~ corresponding to the unrenormalized interaction U. 
We see that for N > AO the values of U PP (N), U hh {N) and 
U p h{N) do coincide and become independent of A^ so for 
the low energy excitations one can define a unique renor- 
malized interaction via U = limjv-^oo U a , a '(N), where 
a, a' = p, h. What is more, because the model with these 
chosen parameters corresponds to the Kondo regime, the 
value of U coincides with the value of irA, clearly seen in 
the inset in figure 1, and as can be deduced from equation 

© for Xc.imp = 0. 



2.1 Results for Symmetric Model 

The estimates of A and U for the symmetric model can be 
checked indirectly from the exact Bethe ansatz results for 
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NRG 


U/irA = 1 


U/irA = 2 




A/A 


U/irA 


R 


A/ A 


U/irA 


R 


A = 2.0 


0.6288 


0.4779 


1.7600 


0.2389 


0.2298 


1.9620 


A = 2.5 


0.6287 


0.4777 


1.7599 


0.2390 


0.2300 


1.9621 


A = 3.0 


0.6288 


0.4780 


1.7601 


0.2390 


0.2299 


1.9619 


A = 3.5 


0.6288 


0.4780 


1.7601 


0.2390 


0.2299 


1.9620 


Exact (BA) 


0.6289 


0.4780 


1.7601 


0.2392 


0.2301 


1.9620 



Table 1. The renormalized parameters for the symmetric Anderson model for U/nA = 1.0 and U/irA — 2.0 calculated using 
the NRG and different values of the discretization parameter A. These results are compared with the corresponding values 
deduced from exact Bethe ansatz results. 



this model |12ll3j . The exact RPT results for the impurity 
spin susceptibility and specific heat coefficient at T = 
for the symmetric model from @ and JSJl are 

1 / U \ 2ir 

Xs ™ = ^a{ 1 + ^a)> 7imp = S' (20) 

By equating these to Xs,imp and 7i mp from the Bethe 
ansatz we can deduce A and U. 

Though the discrete model (|11|) with A > 1 and the 
original model Q with a continuous spectrum have essen- 
tially the same low energy spectrum, i.e. they belong to the 
same universality class, the dependence of the renormal- 
ized parameters on the parameters of the bare model may 
differ. Such differences can occur in calculations where the 
high energy excitations, or high cut-offs, are treated dif- 
ferently. This situation occurs for the iV-fold degenerate 
models (U = oo) where the imposition of the high energy 
cut-off D 1 in the Bethe ansatz calculations for the linear 
dispersion model differs from the band width D for the 
conventional model, but a relation between these cut-offs 
can be found such that the results from the Bethe ansatz 
calculations can be translated into those for the conven- 
tional model |15| . A similar situation applies here. The re- 
sults can be made equivalent by replacing the bandwidth 
D of the discrete model by DAa, with A a is given by 

A A = ±±±^hxA, (21) 

where A a — > 1 in the continuum limit A — > 1. This result, 
which is given in the paper of KWW [2] and also used in 
the paper of Sakai, Shimizu and Kasuya ^7], is derived by 
the requirement that the low energy spectrum of the dis- 
crete and continuum model coincide for U = 0. It means 
that V 2 must be increased by a factor Aa when making a 
comparison with the results of a continuum model with a 
given U/irA. 

The values of the renormalized parameters A, U and 
the xh or Wilson ratio, R = 1 + U/irA |1I2| . deduced 
from these are shown in table 1, where they are compared 
with those deduced from the Bethe ansatz results. The 
agreement with the Bethe ansatz results is remarkably 
good, with errors only of the order of 0.1%, even for values 
of A as large as 3.5. It is important for such accurate 
agreement, particularly for larger values of A that the Aa 
factor is taken into account. 



-0.4 





/ / A/a\ \ 4T k 




/ U/jtA 





I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-3 -2-10 12 3 

U/jtA 



Fig. 2. Renormalized parameters U /irA and A/ A as a func- 
tion of U/irA. The two curves are the results deduced from the 
Bethe ansatz, and the points marked with a cross are those de- 
duced from NRG calculations with A = 2.0. For \U\/wA > 2.0 
the two energy scales merge and asymptotically approach 4Tk 
(dotted line), where Tk is the Kondo temperature. 



In Figure 2 we make a more extensive comparison of 
the results for A and U deduced from the NRG calcula- 
tions and those deduced from the Bethe ansatz for both 
positive and negative values of U. The agreement is ex- 
cellent over the whole range. For U/irA > 2 the energy 
scales merge as the impurity charge becomes localized 
(Kondo regime) and the impurity charge susceptibility 
tends to zero. In this limit U/irA — > 1 and as a con- 
sequence the Wilson ratio R — > 2. The single renormal- 
ized energy scale in this regime is the Kondo tempera- 
ture Tk, and irA = U = 4Tk- In the negative U regime 
such that U/irA < —2 a local bipolaron forms such that 
the spin susceptibility tends to zero, and the charge sus- 
ceptibility is enhanced by a factor of 2. In this limit the 
energy scales again merge such that U/irA — > — 1. This 
is because a Kondo effect develops in the isospin chan- 
nel (doubly occupancy of the impurity site corresponding 
to up-isospin, and zero occupancy to down- isospin), such 
that irA = — U = 4Tk, as the real spin fluctuations are 
suppressed and Xs,imp — * 0. There is a possibility that 
a Kondo effect of this type might be seen in degenerate 
atomic gases with the doubly occupied paired states cor- 
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responding to molecules |16| . The general expression for 
7k that covers both Kondo regimes is 



\U\ 



2\U 



1/2 



-ir\U\/8A+-xA/2\U\ 



(22) 



The analysis of the irrelevant terms about the fixed 
point is not the only way to deduce some of the renormal- 
ized parameters from the NRG calculations. If the NRG 
approach is used also to calculate the dynamics of the im- 
purity model |18I19 | then id and A can be deduced from 
the self-energy and its derivative for u = and substi- 
tuted in equation J5J). Typical values for A calculated in 
this way for the symmetric model (e^ = 0) for U/ttA = 1 
and U/tvA = 2 are 0.6155 and 0.2350, respectively for 
A = 2. They compare well with the Bethe ansatz values 
given in Table 1 . The errors are greater than those deduced 
from the analysis of the fixed point, of the order of 2-3%, 
but still quite small, and are of the order to be expected 
in the calculation of dynamical quantities. Unfortunately 
the evaluation of U from equation © requires a knowledge 
of the four- vertex 7^^(0,0,0,0) which is very difficult to 
calculate directly from an NRG calculation, as it involves 
the Fourier transform of the two-particle Green's function 
with respect to three independent frequency parameters. 
Both A and U have, however, been calculated from equa- 
tions J2J and @ from perturbation theory to third order 
in U for the symmetric model 0], and the results are 



and 



Ail 



U = Uil- 



7T 

T 



7T 2 -9 



u 

ttA 



U 
nA 



(23) 



(24) 



These perturbational results are in good agreement with 
the exact results in the range U/irA < 0.5; at U/ttA = 0.5 
the error in A is less than 1.5% and that for U less than 
4%. 

We were able to define running renormalized parame- 
ters, such as A(N) and id(N) as a function of N, which 
raises the possibility that they can be translated into effec- 
tive parameters appropriate for calculations on an energy 
scale ujn = riDA"^^ 1 ^ 2 or temperature scale, Tn — ujn, 
where r\ is an appropriately chosen constant of order unity. 
We do not however have a unique value for U(N). Never- 
theless, for the particle-hole symmetric case if we take the 
value of Upp(N) (=Uhh(N)) as U(N), and translate this, 
together with A(N) and £d(N), into parameters appropri- 
ate for a temperature scale Tjv, we can generalize the RPT 
expression for the impurity susceptibility in equation Jfjjl 
to finite temperatures, 



Ximp(7 1 ) 



x[° ) p (T)(l + 2^ ) P ( T )) = 



(25) 



where x|m P (T) ^ s tne ^ ree quasiparticle contribution to the 



impurity susceptibility given by 



AO) 

A imp 



(T) = ~\ 



Pimp J— duj 

OLU 




(26) 



Fig. 3. The temperature dependent impurity susceptibility 
Ximp(7)/xim P (0) (stars) versus ln(T/Tk) evaluated from equa- 
tion (I25H for the symmetric model for U/irA — 6.0, compared 
to the results of the Bethe ansatz solution (continuous curve) 
for the s-d model given by Tsvelik and Wiegmann |I2| . 



in units of (g^e) 2 , where f(u) — l/(e"/ T + l), and y5i mp (u;) 
is the free quasiparticle density of states given by equa- 
tion (JJJ. We calculate xl^ p (T) for the symmetric model 
in the Kondo regime for U/ttA = 6.0 at values of Tn, us- 
ing the a renormalized parameters A(Tjsr) and 6^(7^) = 
in evaluating the free quasiparticle density of states for 
this energy scale, with r] = 1.2 as is used in the NRG 
evaluation of spectral densities on a scale ujn (see for ex- 
ample |18I19| ). We then deduce Xi mp (7 1 ) from equation 

(|2*5|) using U(Tn). In figure 3 we compare the results of 
this calculation with the Bethe ansatz results for the s-d 
model given in reference |12|. There is quite a remarkable 
agreement with the exact Bethe ansatz results over this 
temperature range, and the value of Ximp(7 1 ) in the ex- 
treme high temperature range corresponds to that of the 
free bare model, 1/8T. The agreement is much less good 
if we use U p h{N) for U(N) and, as there is not a unique 
prediction for this quantity, one cannot place too much 
reliance on this calculation. However, it does suggest that 
one might be able to define a renormalized perturbation 
theory with running coupling constants (this is possible in 
a magnetic field, with field dependent parameters, as we 
will demonstrate elsewhere). 



2.2 Results for Asymmetric Model 

We now consider how the renormalized parameters vary 
as the bare parameters change from one qualitatively dif- 
ferent regime of the model to another. We start first of all 
with the value td = —~kA = constant, and then increase U 
from U = to U = 2.7ttA. This takes us from what might 
be called the full orbital regime to the Kondo regime. The 
results are plotted in figure 4. We see that id increases at 
first approximately linearly with U until U/irA ~ 1, and 
then slowly increases monotonically through zero at the 
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Fig. 4. The renormalized parameters A/ttA, Id /it A, and 
U/ttA, the Wilson ratio R = 1 + Up~d{Q), and the impurity 
occupation value ni mp are plotted for the asymmetric Ander- 
son model for a fixed impurity level = — n A as a function 
U fix A. Over this range the state of the impurity changes from 
the full orbital to the Kondo regime. 



symmetric point U/ttA ~ 2, remaining very close to zero 
in the Kondo regime at higher values of U. The renor- 
malized resonance width A decreases monotonically over 
the same range, though only slowly at first in the full or- 
bital regime, and approaches zero in the limit U — > oo. 
The quasiparticle interaction U increases at first linearly 
with U, reaching a maximum for U/ttA ~ 1, and then 
decreases so that its energy scale merges with that for A, 
with U = it A = 4Tk, as for the symmetric model dis- 
cussed in the previous section. Over the same parameter 
range n- lmp decreases from an initial value of 1.8 to slightly 
below unity at U — 2.7ttA. The Wilson ratio or \/l ra- 
tio, R = 1 + Uj5d{0), increases from 1 and asymptotically 
approaches 2 in the Kondo regime. For negative values of 
U, Id and U decrease linearly as U decreases and AjA 
approaches unity. In full orbital regime U /ttA < 0.5, in- 
cluding the range with larger negative values of U, the 
Hartree-Fock theory, where id = td + Und/2 and A = A, 
constitutes a reasonably good approximation. 

It is interesting to compare this behaviour with that 
for the spectral density pd{u) of the d-electron Green's 
function, which can be calculated from the NRG results. 
Some of the results for pd(w) over the same parameter 
range are shown in figure 5(i). Initially there is only a 
single resonance which moves to higher energies as U is 
increased. As the Kondo regime is approached this splits 
into a three peaked structure, the central narrow peak at 
the Fermi-level being the many-body Kondo resonance. 
The renormalized Anderson model given by equation Q 
describes only a single resonance, but is valid for the low 
energy behaviour in all the parameter regimes. We see that 
for smaller values of U , and in the negative U regime, id 
tracks at first the lower resonance, increasing monoton- 
ically, and in the Kondo regime tracks the Kondo reso- 
nance. The precise nature of this tracking is made evident 




Fig. 5. Plots of the spectral density pd(u) as a function of a; : 
(i) for values of U fix A with ed = — irA over the same range of 
values as in figure 4, and (ii) for values of ed/irA with U = 2nA 
over the same range as in figure 6. The heights of the maxima 
in each case are normalized to unity so that the position of the 
peak can be compared with the values of id (dotted line) taken 
from figure 4 for plot (i), and from figure 6 for plot (ii). The 
values of U for the curves in (i) can be read off from the left 
hand scale from the coordinates of the corresponding maxima, 
and those for ed in (ii) from the right hand scale. 



in figure 5(i) where the curve for id taken from figure 4 is 
plotted with the spectral density results, where the main 
peak for each Pd{u) has been normalized to unity for the 
comparison; the maxima of these peaks all lie on the curve. 
The height of the atomic- like peaks in this plot, which ap- 
pear at u> ~ ed and u> ~ ed + U for larger U in the Kondo 
regime, become somewhat flattened in the normalization 
as the Kondo peak is so high, making them difficult to see 
clearly in figure 5. 

In figure 6 we give the renormalized parameters for a 
complementary scan. In this case we start first of all with 
ed = —3irA and U — 2ttA, which is in the almost full 
orbital regime as e.d + U < 0, and then increase ed with a 
fixed value for U to a final value ed = I.BttA. In doing so 
we move to a mixed valence regime for e.d + U ~ 0, then 
for ed = —ttA the Kondo regime for the symmetric model 
with U/nA = 2, another mixed valence regime for ed ~ 0, 
and finally an empty orbital regime for ed > ttA. In the 
earlier renormalization group analysis of KWW , these 
parameter regimes are interpreted in terms of essentially 
three different fixed points, (i) the full (empty) orbital (ii) 
the strong coupling and (iii) the mixed valent. In the anal- 
ysis here, however, they are simply different parameter 
regimes of a single renormalized Anderson model. We see 
from figure 6 the characteristic features of full and empty 
orbital regimes, id increases linearly with ed, the Wilson 
ratio R ~ 1, n; mp ~ 2 or 0, and A and U are independent 
energy scales. In contrast in the Kondo regime, id ~ 
and is largely independent of ed L the Wilson ratio R ~ 2, 
fiimp ~ 1, and the energy scales A and U have merged such 
that U = it A. The mixed valence regimes appear more as 
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cross-over regions between these two types of behaviour, 
where the independence of the two energy scales A and U 
emerges. The curve for id from figure 6 is plotted in figure 
5(ii) with the corresponding results for the spectral density 
Pd{^)i with the height of the peaks normalized to unity. 
The renormalized level id is found to track the single peak 
in the full orbital regime as it moves up through the inter- 
mediate valence regime e ( i + U ^ 0, where a second lower 
peak develops near u> ~ e^, and then it tracks the central 
narrow resonance at the Fermi-level in the Kondo regime, 
where a third upper peak develops near lj ~ ea + U . The 
reverse process occurs as €d increases from Cd ~ —irA, as 
the side peaks disappear and a single peak emerges in the 
empty orbital regime. As in figure 5(i) atomic-like side 
peaks at u) ~ e^, and lu ~ td + U, which appear in the 
mixed valence and Kondo regimes, are only just discern- 
able due to the normalization of the height of the central 
resonance. If we concentrate on the mixed valence regime 
for td ~ in figure 6, we see that there is a significant 
upward shift of id- This reflects the effective shift in the 
bare level €d obtained by Haldane [2D, in a poor man's 
scaling treatment in which the virtual charge fluctuations 
were eliminated to focus on the mixed valent regime. For 
the U = oo limit, the thermodynamic behaviour of the 
model in this regime was shown by Haldane to depend 
on the ratio id/ A, where id = td + A/idviirD /2A). For 
the renormalized Anderson model the equivalent ratio is 
id/ A, so id is not to be equated with id but id/ z. In the 
mixed valence regime + U ~ the Haldane shift can be 
seen to be in the opposite direction, to lower energies. 

In the full and empty orbital regimes again the Hartree- 
Fock is a reasonable approximation, particularly when the 
effective level lies away from the Fermi-level, even though 
in this case we have a large value of U. With = — 3ttA, 
U = 2ttA, rid = 1-82, the Hartree-Fock estimate of the po- 
sition of the renormalized level id/nA = (ed+Und/2)/TrA = 
— 1.17 is in good agreeement with the renormalization 
group results, as can be seen in figure 6. The fact that 
the quasiparticle interaction U is comparatively large in 
the full and empty orbital regimes does not imply that 
the interaction effects on the low energy scale in these 
regimes are large. On the contrary they are small, as can 
be seen from the fact that R approaches unity in these 
two regimes. The reason is that the significant term is 
the combination £7/5(0), and in these regimes p(0) — > as 
the renormalized level moves well away the Fermi-level, 
such that U p(0) is quite small. A similar situation applies 
in the case of the U = oo TV-fold degenerate Anderson 
model where slave boson mean field theory is asymptoti- 
cally exact in the limit N — > oo. The mean field theory in 
that case is valid not because U — > as N — > oo, in fact 
U remains finite in the limit, but because p(0) — » so the 
combination (7/5(0) — > as N — > oo |21j . 

Though we have derived a consistent picture in the 
form of a renormalized Anderson model to describe the 
low energy behaviour within both the NRG and RPT ap- 
proaches, the two types of calculations involve quite differ- 
ent ways of realising this renormalization; it is interesting 
to compare and contrast the ways in which they arrive 
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Fig. 6. The renormalized parameters A/nA, id/irA, and 
JJ/ttA, the Wilson ratio R and the impurity occupation value 
fiimp are plotted for the asymmetric Anderson model for a 
fixed value U = 2nA as a function of ed/nA. On increasing 
e<i through this range the impurity state passes from the full 
orbital state, through a mixed valence and Kondo regime, a sec- 
ond mixed valence regime, and to a final empty orbital state. 

at the same results. In the NRG approach the excitations 
at the fixed point are free quasiparticles, and their energy 
spectrum depends solely on the ratio id/ A, and hence 
via the Friedel sum rule, equation (J0J, on the additional 
charge at the impurity site. The leading irrelevant terms 
that cause a scattering of these quasiparticles are of the 
order yl - ^ -1 )/ 2 for large TV [22] and are of two types; 
a one-body or hybridization term which determines the 
width A of the quasiparticle resonance, and a local two- 
body interaction U. This latter term only contributes to 
states in which two or more quasiparticles are excited. 

In the RPT approach the free quasiparticles already 
contain the effective hybridization, and are characterized 
by the two parameters, id, the renormalized impurity level, 
and A, the renormalized width. As in the NRG approach 
the local two-body interaction term U only comes into 
play when two or more quasiparticles are excited. This is 
because the interaction term in the quasiparticle Hamil- 
tonian is normal-ordered with respect to a vacuum corre- 
sponding to the interacting ground state, in contrast to 
the interaction term in the original bare model Q . In the 
more general renormalization prescription of the RPT, 

£(o,o) = o, £'(o,o) = o, f u (o,o,o,o) = tr, (27) 

where E(u>,T) is the renormalized self-energy and the 
renormalized 4-vertex /^(wi, ui2, W3, 0J4) |8I4| . this nor- 
mal ordering is effectively carried out by the counter terms 
that ensure these conditions are satisfied. 

The two approaches also differ in detail for analytic 
calculations for the behaviour of the model at higher tem- 
peratures and higher energy scales. In the NRG approach 
the next to leading order irrelevant terms about the fixed 
point have to be added to the effective Hamiltonian, and 
some of these will correspond to multiple quasiparticle 
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interaction terms. In contrast in the RPT approach, no 
higher order quasiparticle interaction terms have to be 
added; though the counter terms have to be fully taken 
into account order by order in the expansion in U. The 
counterpart of the multiple quasiparticle scattering terms 
are the higher order Feynman diagrams in the pertur- 
bation expansion, which have to be taken into account 
to probe the higher energy scales. Calculations on these 
higher energy scales should describe the 'undressing' of 
the fully renormalized quasiparticles. Calculations to third 
order in U have been carried out for the particle-hole sym- 
metric model 0]. 



3 Model with a coupling to Phonons 

One possible modification of the impurity model is the in- 
clusion of an interaction with phonons. Such a term could 
be important when dealing with an impurity with a par- 
tially filled 4f shell, as there can be a change of the ionic 
volume of the order of 10% when an electron is removed 
or added to the 4f shell due to the adjustment of elec- 
trons in the outer d and s shells. The simplest type of 
interaction to consider is a coupling of the occupation of 
the d or f shells to an Einstein phonon of frequency ujq, 
which is of the form used to study polaronic effects in 
a tight-binding model by Holstein [231 ■ If we add such a 
term to the Hamiltonian for the Anderson model we get 
the Andcrson-Holstein model, which we have studied in 
earlier work [21]. The additional term has the form, 

tf e _ ph = g(tf + b)($2 d Ua - 1) + "O^b, (28) 
a 

where b and W are creation and annihilation operators 
for the phonon modes and a linear coupling has been as- 
sumed with a coupling constant g. Some insight into the 
behaviour of this model can be obtained by performing 
a displaced oscillator transformation, so that the trans- 
formed model takes the form, 

H' =L>" 1 HU = J2 e d,c&dU* + U cS n d ^n dl + 

a 

+ £ V k (e^^dtc ka + e-^^cLA) + 

k,a 
ka 

(29) 

where 

<y = e -4( bt - b )E„«^-i) ; (30) 

and ^ 

Erf.efE = e<i+ — j U cS = U-—. (31) 

LO Ll> 

In the large frequency limit luq — > oo, such that g 2 /u)o 
remains finite, the model is equivalent to the original An- 
derson model with an effective level ed,off and an effective 
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Fig. 7. The values of the renormalized parameters A/ A and 
U/A for the symmetric model with phonons (U — 0) plotted 
as a function of U E {f/irA — — 2g 2 /ujq-kA. The dotted curves 
correspond to the limit wo — * oo (U e s finite) and are the same 
as those for the negative-17 model shown in figure 2. The full 
lines are for the corresponding curves with luq = 0.05. 



local interaction U c g, because in this limit the exponen- 
tial terms that modify the hybridization term are equal 
to 1. When the original Coulomb interaction term [7 = 0, 
the effective interaction is attractive (i7 e ff < 0), so that in 
the localized strong coupling regime local bipolarons form, 
which is equivalent to the formation of a local moment in 
the spin model, and there is a Kondo effect in the charge 
(or equivalently isospin) channel. 

For the large ujq limit, such that we can neglect the 
effect of the exponential terms in the hybridization in the 
Hamiltonian l|29|l . we require ujq to be much greater than 
the conduction bandwidth D [22|, which is unrealistic. For 
the regime wo < O, the physically appropriate one, the 
exponential terms in the phonon operators in (|29|l modify 
the form of the original Anderson model. Nevertheless the 
low energy fixed point, and the leading irrelevant interac- 
tion terms, should correspond to a renormalized Anderson 
model. The explicit inclusion of the phonon terms should 
only be necessary in considering higher order irrelevant 
terms, which contribute only at higher energy scales and 
higher temperatures. The derivation of the renormalized 
parameters in equation (J2J holds provided that the self- 
energy S(u>) is analytic at u> — 0. We can test this con- 
jecture explicitly by again fitting the lower energy levels 
obtained in the NRG calculations to a renormalized An- 
derson model with parameters e<j, A and U. We first of all 
look at the results for the particle-hole symmetric model 
with [7 = 0, where the phonons induce a negative-?/ term, 
and assume particle-hole symmetry to compare our results 
with those obtained in the previous section. In the large ujq 
limit these will be identical to those obtained earlier with 
U being replaced by U e s = U — 2g 2 /uiQ. For comparison we 
consider the results for a value of ujq <c D and a commen- 
surately smaller value of the coupling g, so that it covers 
a similar range of values of U c s ■ The results of such a fit- 
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Fig. 8. The values of the renormalized parameters A/ A and 
17/ A for the symmetric model with phonons, and U/irA = 
2.35, plotted as a function of U c ff/irA = (U — 2g 2 /luq)/itA. 
The dotted curves correspond to the limit too — > oo (U e s finite) 
and are the same as those for the negative-t/ model shown in 
figure 2. The full lines are for the corresponding curves with 
ujq = 0.05. 



ting, using the procedure outlined in the previous section, 
are shown in figure 7 for ujq/D = 0.05, and for compari- 
son also the case corresponding to ujq — > oo. These values 
cover the negative U e s regime as we have U = 0. We see 
that, as the coupling g increases, the values of A decrease 
much more rapidly in the smaller phonon frequency case 
than in the large ojq case due to the exponential terms that 
modify the hybridization in equation H29|) . The values of U 
are also commensurately smaller because in the localized 
or local bipolaron limit, we must have U = —nA, which 
is the case for both sets of results in the strong coupling- 
limit. 

In the results shown in figure 8 we start with a value 
of U/irA = 2.35 for which the impurity charge is almost 
localized and we are in the Kondo regime: as we increase 
g one effect is to decrease the value of U e g, and in the 
large u>o limit this is the sole effect, so that when U e g = 
we are left with an unrenormalized width parameter 
A = A. In the low ojq case the renormalized interaction U 
changes sign at the same point as the bare parameter U e e 
changes sign. This result is somewhat surprising. As the 
effective interaction induced by the phonons is only valid 
on the scale w < wo, one might have expected that the 
on-site interaction U would be significantly renormalized 
on reducing the energy scale to w ~ ujq, say to U, so 
that U would change sign when U — 2g 2 /luq, which would 
be significantly shifted from the condition U c b ~ 0. This 
is not the case which implies that both the direct and 
phonon induced interaction terms are renormalized in a 
similar way. 

At the point where [7 = there are polaronic effects, 
which are absent in the large ojq limit, such that A is 
significantly renormalized. For the values used here the 
width is reduced by factor of approximately 2 when U = 0. 




2g/Q)7tA 



Fig. 9. A plot of the renormalized resonance width A/ A for 
the symmetric model with a phonon coupling g as a function 
of 2g 2 /luq-kA with values of U chosen such that U c « = U — 
2g 2 /loq = 0, for the two cases loo = 0.05 and loo = 0.0125. 



In the negative-t/eff regime the behaviour is as in figure 
7, with bipolaronic localization and U = —ttA, with very 
much reduced values of U and A for the small frequency 
case as compared with the corresponding values for the 
large ujq limit. 

In the corresponding spinless model with a coupling to 
phonons there is no induced attractive on-site interaction, 
but polaronic effects occur due to the retarded effective 
potential 23^ ■ In the model with spins, we can reveal these 
polaronic effects if we choose values of U which cancel the 
attractive interaction that leads to bipolaron formation, 
U e g = U — 2g 2 /loq = 0. We look at the renormalized res- 
onance width for the symmetric model for a range of val- 
ues of the electron-phonon coupling g with appropriately 
chosen values of U. As a consequence the renormalized 
quasiparticle interaction U is negligible. The results for 
the symmetric model are shown in figure 9 for ujq — 0.05 
and luq = 0.0125 (D — 1). It can be seen that the decrease 
in A with increase of g is greater in the smaller phonon 
frequency case. In the limiting case ujq ~~ * °o (2g 2 /wo fi- 
nite) there is no polaronic effect and A/A = 1 in all cases. 
Some insight into the behaviour of the asymmetric model 
(U = 0) with a coupling to phonons can be gained by 
comparing with the negative-J7 c ff Anderson model as it 
corresponds to the case loq — > oo. The asymmetric model 
with negative-[/ can be mapped into a positive-C/ symmet- 
ric model with a finite magnetic field H ~ — under the 
interchange of charge and spin; the one-electron spectral 
density for the isospin up electrons is given by Pd{^) and 
that for the isospin down by pd(—Lu). A finite td favours 
either the doubly occupied state (isospin up) for ed < 0, or 
the empty state for > 0, so that increasing the value of 
| /Jeff | bas the effect of further increasing the polarization. 
The results for the renormalized parameters for the asym- 
metric model with phonons (U = and = — irA) and fi- 
nite wo (= 0.05) which are plotted in figure 10 show a sim- 
ilar trend. As g 2 increases, or equivalently U e s decreases, 
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Fig. 10. The renormalized parameters, A/nA, e^/nA, and 
U /nA, the Wilson ratio R and the impurity occupation value 
fiimp plotted for the asymmetric model ((7 = 0) with a coupling 
to phonons for a fixed value of = — nA as a function of 
Ucs/ttA for lj = 0.05. 

both id and U decrease linearly, and the occupation of the 
impurity level slowly increases and tends to the value 2. 
The interaction effects on the low energy behaviour can 
be seen to be relatively small because, as id moves further 
from the Fermi-level, the quasiparticle density of states at 
Fermi-level p~d(0) decreases and R ~ 1. This trend can be 
clearly seen in the calculated values of the spectral density 
Pd{u) |24| . and id tracks the peak in Pd{^)- The peaks in 
Pd{u) become asymmetric for strong electron-phonon cou- 
pling, as phonon side-bands excitations are induced on the 
low energy side of the peak, giving the appearance of a 
broadened peak. The peak of the renormalized Anderson 
peak, however, is symmetric with a width A slightly re- 
duced from the bare value A, but is quite compatible with 
the spectra seen in Pd{w) for small to. The renormalized 
model has to reproduce the quasiparticle spectrum only 
in the immediate neighbourhood of the Fermi-level which 
in this case lies on the higher energy side of the peak in 

4 Model with a Falikov-Kimball Interaction 

Apart from the interaction U between the electrons in the 
local d- or f-orbitals at the impurity site there will also 
be a two-body interaction between the electrons in these 
orbitals and the conduction electrons. Such a term does 
not appear in the standard Anderson model but a local 
interaction of this type is included in the Falikov-Kimball 
model • It is the sole two-body interaction term in this 
model, which was put forward as an appropriate model for 
investigating valence instabilities in compounds, where the 
occupation of the localized f-orbitals changes significantly 
as a result of pressure or alloying (see reference |27| for a 
recent review). The reason that it is not usually included 
in impurity models, is not because the interaction is very 
small, but because it is thought not to play an essential 
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Fig. 11. The renormalized parameters A/itA, Id/nA, and 
JJ/ttA, the Wilson ratio R and the impurity occupation value 
fiimp are plotted for the Anderson model with an additional 
Falikov-Kimball interaction Uik- The values of U and Ufk are 
fixed, such that f/fk = U/2 = irA, and ed/irA is varied over the 
range indicated. Also shown as dotted curves are the renormal- 
ized values for the Anderson model with (7fk = 0, correspond- 
ing to those shown in figure 6, but with a shift of such that 
e<i —>■ £d — f/fk- There is a remarkable agreement with these 
displaced results over this range. 

role in understanding the impurity behaviour. This means 
that its effects can, for the most part, be absorbed as a 
rcnormalization of terms that are already within the stan- 
dard models. We can put this hypothesis to the test by 
including such a term in the Hamiltonian and then con- 
sider what effect it has on the renormalized parameters 
of the low energy fixed point. In the discrete linear chain 
version of the Anderson model, equation Qllfl. used in the 
NRG calculations such a term can be included as an inter- 
action between the impurity and the electron occupation 
at the first site on the conduction electron chain, 

ct : ct' 

where (7fk denotes the matrix element of the interaction. 
We add this term to the Hamiltonian in equation (|ll|l and 
then examine the nature of the low energy fixed point. We 
take U = 2irA, as earlier, with Ufk = U/2 = irA, and take 
a range of values of td to run from the full orbital, through 
the Kondo, to the empty orbital regime. We can then com- 
pare the results with those obtained from a similar study 
in section 3. 

The results for the renormalized parameters in this 
study are shown in figure 11. The low temperature be- 
haviour of the model, x s ,imp, Xc.imp and 7i m p, can be de- 
duced by substituting the renormalized parameters into 
equations © and © . There is a clear correspondence with 
results derived earlier for the Anderson model without this 
additional term, shown in figure 6. To examine the cor- 
respondence in more detail we have replotted the values 
from figure 6 in the same figure, using dotted lines, and 
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have displaced the curves by a constant shift of irA = Ufa, 
such that the original bare level td from figure 6 is replaced 
by an effective bare level Id = ej, + Ufa. The agreement 
between the results with the additional Falikov-Kimball 
term, Ufa — irA, and those of the model with Ufa = 
and an effective bare level Zd are quite remarkable. In the 
strong correlation regime the differences are less than 5%, 
and overall less than 10%. As the nature of the low en- 
ergy fixed point has not been changed by the inclusion of 
the Falikov-Kimball term, it is not surprising that the low 
energy behaviour can still be described by a renormalized 
Anderson model. What is somewhat unexpected in this 
case is that these renormalized parameters correspond to 
a constant shift of the bare level over the whole parameter 
range, and that no adjustment of the bare hybridization 
term is required. 

To examine whether this shift varies or not with U , we 
take €d = —2nA with Ufa — it A so that + Ufa = —nA, 
and vary U over the same range as in the results for the 
pure Anderson model shown in figure 4. If the effective 
bare level id is independent of U the results should co- 
incide with those shown in figure 4. The results of these 
calculations are shown in figure 12. Again we see that the 
results are close to those of the Anderson model without 
this term, with the same shift in the bare td- The difference 
in values are less than 4% for U/irA > 1, and less than 
12% overall. In these results, and those shown in figure 
11, the differences in the values of ni mp are even smaller, 
less than 1%, over the whole parameter range. 

We conclude that the low temperature behaviour of 
the model with a Falikov-Kimball interaction Ufa in this 
parameter regime is very similar to that of an Anderson 
model without this term, but with a shift of e<j to Id- The 
shift is approximately independent of td and U, and linear 
for Ufa in the range < Ufa/nA < 2.0. In this range we 
found the deviations between the results of the effective 
model with Id and the model with the Falikov-Kimball 
term increase linearly with the value of Ufa. 

In the case of a vanishing hybridization the nature 
of the fixed point Falikov-Kimball term will be different 
as it will correspond to an X-ray type of problem, with 
the localized impurity electron corresponding to a core- 
hole, and the self-energy at ui — will be singular |28) . 
This, however, is a special case, and in general we can 
expect the self-energy to be regular at u> = and the 
renormalized Anderson model to provide the appropriate 
low energy model (see also the discussion in Sec. V. B of 
reference 

5 Conclusions 

We have shown that the two renormalization approaches 
to the calculation of the low energy excitations of the An- 
derson model lead to a consistent picture of the low en- 
ergy behaviour of the model as described by a renormal- 
ized version of the same impurity model. It is interesting 
to make a comparison of this situation with field theoretic 
calculations for the </> 4 model or quantum electrodynamics 
(QED). In the perturbation calculations for these models 
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Fig. 12. The renormalized parameters A/irA, e^/nA, and 
U/tvA, the Wilson ratio R and the impurity occupation value 
«i mp for the Anderson model with an additional Falikov- 
Kimball interaction Ufk, using the same value as for figure 11. 
Here, td = — 2irA, and U/ttA is varied over the range indicated. 
These results are compared with the corresponding values for 
the Anderson model with Ufk ~ (dotted curves) shown in 
figure 4, for which = —nA. There is again a close agreement 
with the results of the model without the additional term but 
with a displaced ta- 

infinities arise due to the lack of a high energy cut-off, but 
the RPT, with the introduction of counter-terms, allows 
one to by-pass these infinities and perform calculations 
in terms of a renormalized version of the same model, 
with particle masses and interactions taken from experi- 
ment. The RPT calculations for the Anderson model pro- 
ceed along similar lines, and the perturbation theory be- 
comes one in terms of the fully dressed quasiparticles and 
their interactions, in 1-1 correspondence with the origi- 
nal model. Here, however, the renormalizations, though 
in some situations very large, are not infinite as the con- 
duction band edges provide a high energy cut-off, and the 
renormalized parameters, id, A, the analogue of the renor- 
malized masses can be calculated explicitly from equation 
(01 once the self-energy S(u>) is known. 

The Anderson model, as a model for impurity systems, 
neglects many higher energy scale interactions. In includ- 
ing an electron-phonon coupling and a Falikov-Kimball 
term we have included some of the possible higher energy 
interaction terms. Explicit calculations have shown that 
these terms do not change any essential features of the low 
energy behaviour, which can still be described in terms of 
a three-parameter renormalized Anderson model, though 
the renormalized parameters are modified. It is reason- 
able to assume that if other types of interaction terms 
are included, such as two-body hybridization terms aris- 
ing from off-diagonal elements of the Coulomb interaction 
of the impurity d-electron with the conduction electrons, 
that the same would hold. The fact that we can describe 
the various forms of low energy behaviour of this class of 
non-degenerate impurity models, within the framework of 
a single renormalized model, gives a unifying perspective. 
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It also simplifies the interpretation of experiments on the 
low temperature behaviour of impurity systems if the re- 
sults can be analysed in terms of a single three-parameter 
model. 

The non-degenerate Anderson model is not sufficient 
to describe many impurity systems where the impurity ion 
has localized incomplete d or f-shells. There are general- 
izations of the Anderson model that take account of the 
degeneracy of these configurations, such as the U = oo N- 
fold degenerate model for 4f ionic systems, and the n-fold 
denenerate models with Hund's rule coupling for d-state 
ions, and renormalized versions of these models also de- 
scribe the low energy regime. Explicit calculations of the 
renormalized parameters have only been carried out in the 
Kondo regime, with renormalized parameters, such as the 
renormalized Hund's rule coupling expressed in terms of 
Tk |14l21j . Models that have a non- Fermi-liquid low en- 
ergy fixed point, however, are outside the scope of the 
present discussion. 

Though the extra interactions, such as the coupling 
to phonons or the Falikov-Kimball term, need not be in- 
cluded explicitly to describe the very low energy behaviour, 
the question arises as to whether they are needed to de- 
scribe the behaviour on higher energy scales, or whether 
can they be simply taken into account via a re-parameterizat 
of the bare model to Ed, A and U |21I29| . There is no sim- 
ple answer to this question. There are situations where 
this is formally possible, such as in the model with cou- 
pling to phonons in the o>o — > oo limit, where the phonons 
can be effectively eliminated to give a renormalized bare 
model, as given in equation H31(l . However, for u> > loq the 
real phonon excitations, and not just virtual ones, have 
to be taken into account explicitly. Even the interaction 
term U of the bare Anderson model is a renormalized 
one because its experimental value does not correspond 
to the Coulomb matrix elements for electrons in localized 
d-orbitals, but takes into account the many-body relax- 
ation and screening effects of the other electronic shells, 
when a d or f-electron is removed; so the model is already 
a renormalized or effective model. This is true for most 
bare models: which interactions have to be taken into ac- 
count explicitly is dependent on the energy scales and the 
type of experiment which is being described. 



single particle excitations for U = so that the Green's 
function at the impurity site (i = — 1) is given by 



(33) 



where <7oo(w) is the Green's function at site i = for the 
system described by the Hamiltonian Hq n . This Green's 
function is given in turn by 



SooM 



f 



(34) 



and <7n(o;) is the Green's function at site i — 1 for the 
system described by the Hamiltonian H® N . This proce- 
dure can be extended to express goo{E) in the form of a 
continued fraction. 

The one-particle excitations E are given by the poles 
ofljSSJ), 

E e d A^/ 2 - V 2 A^g 0Q (E) = 0, (35) 
where goo(E) can be expressed in the form, 



ion 



9oo 



l<M0)l 2 



(36) 



where <fii(0) and Ei(N) are the eigenvectors and eigenval- 
ues of Hq n . 

If Ep(N) and E®(N) are the lowest particle and hole 
excitations from the ground state of the Hamiltonian H l ^ 1 N 
then, 



V 2 



-E° h (N)A- 



-(JV-l)/2 



V 2 



V 2 



6rf 
V 2 



A^/ 2 g 00 (E°(N)), (37) 



A^/ 2 g m (-El{N)). (38) 



One of us (ACH) wishes to thank the EPSRC (Grant GR/S18571 
for financial support, and W. Roller for helpful advice. AO 
wishes to acknowledge support by the Grant-in- Aid for Scien- 
tific Research from JSPS. We are also grateful to Y. Shimizu 
and O. Sakai for helpful discussions, and particularly in relation 
to the derivation of the correction factor Aa for the discretized 
model. 



A Appendix 

Let the linear chain Hamiltonian of the non-interacting 
system plus impurity be denoted by 1 N , and that of 
the rest of the chain without the impurity starting at site 
i by Hf N , where i = 0, 1, 2 . . .. We are considering free 



If we replace E®, h <. by the corresponding values for the 
interacting system U ^ 0, E p (h\, then they can be used 
to define iV-dependent quantities, e<i{N), and A(N) = 
nV(N) 2 /2, 

I^ \e p {n)a-^/ 2 _ kUN) = A{N _ iy 

2A(N) 2A(N) 



>,/2 9oo(E p (N)), 
(39) 



2A(N) 2A(N) y V " 

. . (40) 

The renormalized parameters q and A corresponding to 

the low energy fixed point are given by = limjv— »oo ^d(N) 

and A — linijv^oo A(N), and are determined by the two 

equations, 

?H= lim ^- 1 )/ 2 .g o(±^)(^)), (41) 

2A N~*oo 
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and 

_ lim A( N -^(g 00 (E p (N))-g 00 (-E h (N))) 
2 A JV-°° E p (N)+E h (N) ' 1 ' 

The results in equations (|41() and (|42|l are applicable 
also systems with non-symmetric or non-constant densi- 
ties of conduction states, with the appropriate form for 
goo(u>), which can be evaluated using either equation i|36|l 
or 1231). 
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